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ABSTRACT 

A  computer  model  simulating  performance  of  passive  optical  sensors  for  imaging 
underwater  objects  from  above-water  paltforms  is  described.  The  model  is  based  on 
the  models  of  sun  and  sky  spectral  radiance,  wavy  sea  surface,  light  propagation  in 
seawater  and  is  able  to  render  images  of  underwater  objects  as  seen  through  the  air- 
water  interface. 
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Executive  Summary 

A  computer  model  has  been  developed  to  simulate  tire  performance  of  passive 
optical  imaging  sensors,  such  as  a  human  eye,  a  conventional  digital  camera  or  a 
hyperspectral  sensor.  The  model  is  based  on  the  ray  tracing  technique  and  it 
implements  several  basic  models  describing  various  elements  of  image  formation. 
These  models  are  an  analytical  model  of  sun  and  sky  radiance,  a  model  of  wind  driven 
waves  on  the  sea  surface  based  on  the  directional  wave  spectrum,  a  sun  and  sky 
reflection  model  by  capillary  waves  on  the  sea  surface,  and  a  bio-optical  model  of 
seawater.  Several  empirical  and  semi-empirical  relations  are  used  for  describing  such 
elements  as  radiance  and  irradiance  transfer  through  a  wavy  sea  surface,  image  blur 
due  to  light  scattering  in  ocean  water  etc.  An  analytical  expression  for  radiance 
scattered  into  a  ray  path  has  been  developed  in  the  present  work. 

All  these  models  are  based  on  approximations  of  various  degrees  of  accuracy. 
Some  of  them  were  verified  in  the  present  work  against  other  models  and  ground  truth 
measurements.  The  model  produces  realistic  images  of  sea  surface  and  underwater 
targets  but  further  verification  of  the  simulation  model  and  its  components  is  required 
and  planned  in  the  near  future.  The  work  will  include  comparison  of  simulations  with 
hyperspectral  data  and  ground  truth  measurements.  Further  improvements  to  the 
model  elements  will  be  made.  It  is  planned  to  expand  the  model  to  incorporate 
simulation  of  the  performance  of  active  imaging  systems  such  as  Laser  Range  Gated 
(LRG)  cameras  and  the  Streak  Tube  Imaging  Lidar  (STIL).  Also,  atmospheric  effects  on 
light  propagation  near  the  sea  surface  will  be  included  in  the  model. 

The  verified  model  may  be  used  for  estimating  performance  of  optical  sensors  in 
various  scenarios  of  underwater  object  detection.  Output  of  the  model  may  provide 
input  for  such  areas  of  operations  research  as  mine  hunting  and  amphibious 
operations. 
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1.  Introduction 


Maritime  Operations  Division  of  DSTO  is  currently  assessing  various  electro- 
optical  techniques,  both  passive  and  active,  for  underwater  object  detection.  Although 
the  final  word  in  such  an  assessment  belongs  to  experiment  the  initial  estimation  of  a 
system  performance  can  be  made  via  simulation  based  on  a  physically  correct  model  of 
image  transfer  through  water  and  the  air-sea  interface.  Another  advantage  of  such  a 
model  is  that,  once  verified,  it  can  be  used  for  estimation  of  sensor  performance  under 
various  operational  conditions  where  real  experiment  may  be  expensive. 

A  model  is  currently  being  developed  for  simulating  images  of  underwater  objects 
as  seen  by  various  types  of  optical  sensors,  both  passive  and  active.  The  simulation 
model  is  based  on  physically  accurate  models  of  sun  and  sky  radiance,  wind  generated 
waves  on  the  sea  surface,  and  the  optical  properties  of  ocean  water.  The  model  can 
generate  images  of  a  three-dimensional  object  constructed  from  simple  geometrical 
objects  such  as  cylinders,  spheres,  cones  etc.  It  can  also  take  as  input  any  digital  image 
and  simulate  its  degradation  due  to  scattering  in  water  and  through  the  wavy  sea 
surface.  This  report  describes  the  model  implementation  for  passive  imaging  sensors 
and  presents  the  results  of  simulations. 

The  geometry  and  main  features  of  the  passive  sensor  model  are  shown  in 
Figure  1.  The  simulation  of  an  image  is  based  on  the  ray  tracing  technique.  Rays  are 
traced  from  each  pixel  in  the  image  plane  and  radiance  towards  the  receiver  is 
calculated  in  the  ray  path.  The  direction  of  rays  in  water  is  calculated  from  Snell's  law 
applied  to  a  wavy  sea  surface.  Since  the  model  simulates  a  passive  sensor  an  accurate 
model  of  natural  illumination  plays  an  important  role.  The  path  radiance  in  an 
individual  field  of  view  is  the  sum  of  sun  and  sky  radiance  reflected  from  the  sea 
surface,  backscattered  from  water,  reflected  from  the  bottom,  and  reflected  from  a 
target  if  the  ray  crosses  the  target. 
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Image 


Figure  1.  Passive  sensor  model  features.  The  names  of  the  authors  indicate  the  models  used. 

2.  Simulation  Software  Design 

All  model  components  have  been  programmed  and  implemented  in  the  software 
simulating  passive  images  of  underwater  objects.  The  skeleton  of  this  software  is  the 
ray-tracing  loop.  All  rays  are  independent  of  each  other.  They  are  generated  one  by  one 
starting  from  the  image  plane  of  a  sensor.  Then  their  path  is  followed  through 
atmosphere,  sea  surface  and  water,  and  radiance  in  the  ray  path  is  calculated  at  various 
parts  of  the  ray  using  the  corresponding  models. 
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Figure  2.  The  simulation  software  flow  chart. 
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The  simulation  software  is  written  in  C++  in  platform  independent  way.  It  can  be 
run  on  a  personal  computer  or  on  a  UNIX  machine.  The  input  parameters  are  provided 
as  a  text  file  that  can  be  generated  on  a  PC  using  a  Graphical  User  Interface  (GUI) 
developed  specifically  for  this  application. 

The  simulation  software  flow  chart  is  presented  in  Figure  2.  At  the  beginning  of 
the  ray-tracing  loop  a  ray  is  generated  from  a  pixel  in  the  image  plane  using  input 
parameters  for  the  sensor  viewing  geometry.  Then  the  program  determines  whether 
this  ray  is  looking  above  horizon  and,  if  this  is  true,  calculates  radiance  into  the  ray 
from  the  direction  of  sky  pointed  by  the  ray  using  the  sun  and  sky  radiance  model. 
This  radiance  is  the  only  contribution  to  the  total  radiance  field  at  the  sensor  and  we  go 
to  the  next  pixel  to  generate  a  new  ray. 

If  the  ray  direction  is  below  horizon,  we  determine  next  whether  it  crosses  the 
above  water  part  of  the  object,  if  there  is  any.  If  this  is  the  case,  the  program  calculates 
radiance  reflected  from  the  above  water  part  of  the  target  using  the  sun  and  sky  model 
and  the  target  shape  and  reflectivity  parameters.  Again,  this  is  the  end  of  this  ray  and 
we  go  on  to  the  next  pixel. 

If  the  target  is  completely  underwater  or  this  specific  ray  does  not  intersect  the 
target  surface,  we  need  to  find  the  ray  intersection  with  the  sea  surface  using  the  sea 
surface  realisation.  From  the  latter  we  also  find  the  wave  facet  slopes  and  the  surface 
elevation  at  the  point  of  intersection.  Variances  of  capillary  wave  slope  distribution  are 
calculated  from  the  sea  surface  model.  This  information  is  used  to  calculate  the  sun  and 
sky  radiance  reflected  from  the  sea  surface.  This  radiance  is  added  up  to  the  total 
radiance  at  the  sensor  and  the  ray  is  followed  further  through  the  sea  surface.  The 
wave  slopes  calculated  from  the  sea  surface  realisation  are  used  to  calculate  the 
direction  of  the  refracted  ray. 

The  program  then  determines  whether  the  refracted  ray  crosses  the  underwater 
object.  If  it  does,  the  radiance  reflected  from  the  target  surface  is  calculated  using  the 
model  for  the  underwater  light  field.  The  reflected  radiance  is  then  transferred  through 
the  water  taking  into  account  the  model  for  light  attenuation  by  water.  It  is  then 
transmitted  through  the  water  surface  and  added  up  to  the  total  radiance  at  the 
receiver. 

If  the  ray  does  not  cross  the  target,  it  is  followed  to  the  bottom.  The  radiance 
reflected  from  the  bottom  is  calculated  and  treated  in  similar  way  to  that  reflected  from 
an  underwater  object.  It  is  added  up  to  the  total  radiance  at  the  sensor,  and  the 
program  goes  on  to  the  next  pixel. 

Once  the  loop  through  all  the  pixels  is  finished,  the  total  radiance  at  the  receiver  is 
convolved  with  sensor  spectral  sensitivity  to  produce  an  image. 
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3.  Sun  and  Sky  Radiance  Model 


To  simulate  the  image  of  an  underwater  object  we  need  a  reasonably  accurate 
model  of  the  light  field  illuminating  the  object.  For  passive  sensors  such  illumination 
originates  from  direct  sun  light  and  light  scattered  by  the  atmosphere.  Therefore,  we 
need  a  model  of  spectral  sun  and  sky  radiance  at  sea  level.  MODTRAN  is  widely  used 
in  the  scientific  community  for  such  calculations  and  probably  represents  the  most 
accurate  and  universal  model  of  sun  and  sky  radiance.  However,  it  is  time  consuming 
and  it  is  worth  having  a  simpler  model,  which  may  describe  a  limited  range  of 
atmospheric  conditions  but  is  sufficiently  accurate  and  robust. 

One  such  simple  model  is  described  by  Zibordi  and  Voss  (1989).  The  model  is 
based  on  an  approximate  solution  of  the  radiative  transfer  equation  for  a  plane  parallel, 
homogeneous,  cloudless  atmosphere  bounded  by  a  lambertian  surface.  Zibordi  and 
Voss  (1989)  compared  the  model  with  measurements  performed  for  different 
wavelengths  and  solar  zenith  angles  and  found  that  agreement  in  most  cases  is  within 
±10%.  A  detailed  description  of  the  sky  radiance  model  can  be  found  in  Appendix  A, 
which  mainly  reproduces  the  description  from  the  paper  of  Zibordi  and  Voss  (1989) 
and  is  given  here  for  completeness. 

We  compared  this  model  with  MODTRAN  calculations.  A  rigorous  comparison 
would  require  many  time  consuming  MODTRAN  computations  and  we  leave  such  an 
exercise  for  future  work.  Here  we  make  comparison  at  selected  parameters.  The 
following  figures  show  direct  and  diffuse  solar  irradiance  at  a  horizontal  surface  and 
sky  radiance  at  several  viewing  angles.  For  all  computations  the  sun  zenith  angle  was 
45  degrees.  The  Zibordi  and  Voss  model  includes  the  following  parameters:  Angstrom 
parameters  for  aerosol  optical  thickness  a0  =  0.05,  /?  =  0.8,  the  asymmetry  parameters 
for  the  two  term  Henyey-Greenstein  (TTHG)  function  a=0.962,  gi=0.713,  and  £2=0.759, 
ozone  concentration  C0j  =  0.33  cm  and  water  vapour  concentration  Cw  =  2.5  cm.  It 

was  found  by  Zibordi  and  Voss  that  these  parameters  most  accurately  describe  their 
experimental  results.  The  field  experiment  was  conducted  on  fine  day  of  15  May  on 
California  coast.  In  the  MODTRAN  calculations  we  used  default  parameters  for  the 
midlatitude  summer  model  with  navy  maritime  aerosol  model  for  open  ocean.  From 
our  point  of  view  this  represents  the  closest  to  the  conditions  of  experiment  default 
model.  Three  sets  of  calculations  were  performed  with  the  single  scattering  option,  the 
two-flow  multiple  scattering  (MS)  model,  and  the  DISORT  model  for  multiple 
scattering  with  sixteen  streams. 
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Figure  3.  Direct  irradiance  at  horizontal  surface.  Sun  zenith  angle  is  45  degrees. 


Figure  3  compares  direct  irradiance  at  the  horizontal  surface  using  the  two 
approaches.  Reasonably  good  agreement  can  be  seen  in  the  visual  part  of  the  spectrum 
with  most  disagreement  in  the  near  infrared  part  of  the  spectrum,  which  is  not  of  much 
significance  for  underwater  studies. 
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Figure  4  shows  relatively  good  agreement  between  MODTRAN  and  the  Zibordi 
and  Voss  model  for  diffuse  irradiance  at  horizontal  surface  with  the  latter  model  giving 
higher  values  in  blue  part  of  the  spectrum. 

Figures  5  to  7  show  sky  radiance  at  various  viewing  angles.  The  main  conclusion 
from  these  figures  is  that  the  single  scattering  approximation  in  MODTRAN  gives 
quite  poor  results.  Even  the  default  two-stream  multiple  scattering  approximation  is 
not  as  good  as  the  DISORT  option  with  16  streams.  Thus,  if  we  need  to  take  into 
account  sky  radiance  very  accurately,  we  will  need  to  use  the  DISORT  option  in 
MODTRAN  with  at  least  16  streams,  which  will  increase  calculation  time  dramatically. 
The  computation  time  of  one  viewing  direction  for  a  Pentium  II  500MHz  computer  was 
1  min  45  sec  for  single  scattering  mode,  2  min  33  sec  for  the  default  two-stream 
multiple  scattering  option,  and  about  90  min  for  the  DISORT  option  with  16  streams. 
To  represent  radiance  all  over  the  sky  accurately  we  need  about  200  viewing  angles. 
This  equates  to  about  300  hours  of  MODTRAN  computations  for  one  atmospheric 
scenario.  By  atmospheric  scenario  we  mean  a  single  value  of  sun  elevation  angle  above 
horizon  for  a  given  cloud  cover,  aerosol  model  etc.  Precomputed  sky  radiance  can  be 
easily  rotated  by  any  azimuthal  angle,  so  the  sun  azimuth  could  be  set  to  zero  during 
preliminary  computation  of  the  sky  radiance.  We  should  note  here,  however,  that  if  the 
cloud  cover  is  not  isotropic,  it  will  be  fixed  with  respect  to  the  sun  position,  so  sun  and 
sky  can  be  rotated  as  a  whole  around  the  sensor. 


Figure  5.  Sky  radiance.  Sun  zenith  angle  is  45  degrees,  vieioing  zenith  angle  is  5  degrees,  and 
viewing  azimuth  is  180  degrees  from  the  direction  to  the  sun. 
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Figure  6.  Sky  radiance.  Sun  zenith  angle  is  45  degrees,  viewing  zenith  angle  is  45  degrees,  and 
viewing  azimuth  is  180  degrees  from  the  direction  to  the  sun. 


Wavelength  (nm) 


Figure  7.  Sky  radiance.  Sun  zenith  angle  is  45  degrees,  viewing  zenith  angle  is  45  degrees, 
viewing  azimuth  is  7.2  degrees  from  the  direction  to  the  sun. 
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4.  Sea  Surface  Model 

The  sea  surface  is  modelled  in  the  following  form  (Massel,  1996): 

M  N 

=  c°s(^c°s0„  +  fcjsin0„ -wmt  +  q>m„)  (4.1) 

m= 1  n=\ 

Here  amn  is  the  amphtude  of  individual  waves,  km  is  the  wavenumber,  (Om  is  the  wave 
frequency,  and  cpmn  is  the  phase  of  individual  waves.  M  and  N  are  the  number  of 
intervals  in  the  wavenumber  and  azimuth  0 ,  respectively. 

Coefficients  in  equation  (4.1)  are  calculated  from  a  wave  spectrum  using  the 
following  procedure.  First  we  determine  the  range  of  wavelengths  in  the  random 
realisation  of  the  sea  surface.  The  minimum  wavelength  is  set  to  be  twice  the  minimum 
pixel  size  at  sea  surface.  The  maximum  wavelength  is  twice  the  maximum  size  of  the 
whole  field  of  view  footprint  at  sea  surface.  The  optical  properties  of  the  sea  surface 
due  to  shorter  waves  are  modelled  statistically. 

The  wavenumber  interval  [fcmin,  kmax],  determined  by  the  maximum  and  minimum 
wavelength,  respectively,  is  divided  into  M  intervals  with  middle  points  km  and 
interval  sizes  A km.  The  wavenumber  km  is  randomly  chosen  from  each  interval  as 
k  =k  +  (r-  0.5)AJL, ,  where  r  is  a  random  number,  0  <  r  <  1.  The  direction  of  a  mode 

propagation  is  chosen  randomly  as  6n  =6n+  (r  -0.5)A6n ,  l<n<  N  .  The  phase  of  a 
mode  is  a  random  number  in  the  interval  [0,  2k],  (pm„  =  2n r . 

The  amplitude  of  a  mode  amn  is  calculated  from  a  wavenumber  spectrum  ^{k^O  ) 
(Phillips,  1977): 

amn=^(km,0n)knAkmMn  (4.2) 


The  frequency  of  the  mode  is  determined  from  the  dispersion  equation  (LeBlond 
and  Mysak,  1978): 


CO2  = 


gk  + 


ok 


3  y 


tanh  kH , 


(4.3) 


) 


where  g  is  the  specific  gravity,  o  is  the  water  surface  tension,  p  is  the  water  density,  H 
is  the  water  depth.  For  gravity  waves  in  deep  water,  equation  (4.3)  can  be 
approximated  as  follows: 
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CO2  -gk 


(4.4) 


This  approximation  is  accurate  for  most  of  the  applications  of  imaging  through  the 
water  surface  that  we  are  concerned  with  in  the  present  work. 

Any  wave-number  spectrum  can  be  used  to  generate  a  sea  surface  realisation.  In 
our  simulations  we  used  the  directional  wave-number  spectrum  of  Apel  (1994) 
(Appendix  B). 

5.  Radiance  Reflected  from  Sea  Surface 


Reflection  of  sun  and  sky  radiance  from  a  wavy  sea  surface  is  modelled  following 
the  work  of  Zeisse  (1995).  In  general  terms  the  reflected  radiance  can  be  expressed  by 
the  following  equation: 


Lr — j*  r  t»<7r  /  2  p(o>)L,  (e„  <P,  )q(e, ,  <p, ,  f , ,  c , ,  , 

**  w  Ur=const 


(5.1) 


where  pico)  is  the  Fresnel  reflection  at  the  angle  of  incidence  CO,  9r, cpr  are  the  zenith 
and  azimuth  angles  to  the  receiver,  6X ,  (p  t  are  the  zenith  and  azimuth  angles  of  ray 
reflected  from  a  wave  facet  towards  sky,  Ls  (ds  ,(p  j  is  the  radiance  from  sky  or  sun, 
Cx,£y  are  the  slopes  of  a  wave  facet,  q{Or,g>r,^x,^y,W)  is  the  interaction  probability 

density  (Plass  et  al.,  1975),  and  W  is  the  wind  speed.  The  latter  can  be  expressed  via  the 
occurence  probability  density  and  calculated  using  slope  variances  obtained  from  a 
surface  wave  spectrum  for  waves  shorter  than  the  pixel  footprint  size  at  the  sea  surface. 
The  integral  conditions  indicate  that  the  angles  of  incidence  less  than  90  degrees  are 
considered  and  the  direction  to  the  receiver  is  kept  constant  during  integration  over 
possible  wave  slopes  at  the  sea  surface.  In  the  case  of  reflection  of  the  sun,  equation 
(5.1)  can  be  significantly  simplified  due  to  the  small  value  of  the  solid  angle  subtended 
by  the  solar  disk  (see  Zeisse,  1995  for  details). 

6.  Radiance  from  below  sea  surface 

The  calculation  of  radiance  scattered  from  water  in  the  direction  of  the  receiver 
along  a  ray  path  is  illustrated  in  Figure  8.  In  the  computation  of  radiance  we  make  the 
following  assumptions: 

-  we  assume  the  single  scattering  approximation; 

-  the  optical  properties  of  seawater,  both  inherent  and  apparent,  are  uniform  with 
depth; 
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-  the  volume  scattering  function  of  water  is  constant  in  the  direction  of 
backscattering; 

-  the  mean  cosine  of  the  downwelling  light  stream  is  constant  with  depth. 

These  assumptions  simplify  the  analysis  significantly  and  are  sufficiently  accurate 
for  initial  estimation  of  system  performance.  Further  improvement  of  the  model  is  left 
for  future  work.  This  may  require  a  numerical  solution  for  the  light  field  in  water. 


Figure  8.  Calculation  of  radiance  into  individual  pixel's  field  of  view  (FOV). 


Intensity  of  light  scattered  in  the  direction  of  the  receiver  by  the  volume  dv  =  dzds 
is  by  definition  (Jerlov,  1976): 


dl  =  p(a)dv 

VAZ) 


(6.1) 
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where  (5(a)  is  the  volume  scattering  function  of  seawater,  Ed(z)  is  the  downwelling 

irradiance  at  the  depth  2,  jld(z)=  cos  6  d  is  the  mean  cosine  of  the  downwelling  light 

stream  at  depth  z,  a  is  the  angle  between  the  mean  direction  of  downwelling  light  and 
a  ray,  such  that: 

cos  a  =  cos0d  cos  9  +  sin  6d  sin  6?  coscp,  (6.2) 

where  (p  is  the  azimuth  angle  between  the  direction  to  receiver  and  downwelling  light 
field. 

After  some  transformations  we  find  that  radiance  scattered  from  the  volume 
element  into  an  individual  pixel's  field  of  view  (FOV)  just  below  the  sea  surface  is 


rfi>(o-)=AMJ3(a)exp 

Kd+-^~\ 

fij  cos  0 

_ 

cos  6  J 

(6.3) 


Here  we  assume  that  the  radiance  attenuation  coefficient  KL ,  the  mean  cosine  of 
the  downwelling  light  jLLd  and  the  attenuation  coefficient  of  the  downwelling 
irradiance  Kd  is  constant  with  depth,  so  the  downwelling  irradiance  at  depth  z  is 

£</(*)  =£</(0")ex  v(~Kdz),  (6.4) 

where  Ed  ((T )  is  the  downwelling  sun  and  sky  irradiance  just  below  the  sea  surface. 

The  radiance  reflected  from  water  along  the  whole  ray  path  is 


i„(o-)=}<(o-)fe  (6.5) 

0 

We  assume  here  that  the  volume  scattering  function  is  constant  in  the  direction  of 

b  n 

backscattering  ,  i.e.  for  nil  <a  <K  (3(a)  =  — ,  where  bb=2n  [ (5(a)s'mada  is  the 

2  71  ^ 

nil 

backscattering  coefficient.  This  is  a  normal  assumption  in  an  analytical  approach  to 
radiative  transfer  in  seawater  based  on  a  characteristic,  highly  forward  peaked,  volume 
scattering  function  of  seawater  (Zege,  1991).  Assuming  also  that  jld  is  constant  with 

depth  we  have  the  following  equation  for  the  radiance  just  below  the  sea  surface  in  an 
individual  pixel's  FOV 
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L  (0-)  =  ^fck _ 

"v  '  2 K{Xd  KdcosQ  +  K 

The  radiance  just  above  the  sea  surface  is  (Jerlov,  1976) 

l  (0-)=Z»(°  \~P^\  (6.7) 

n 

where  n  is  the  refractive  index  of  sea  water  relative  to  air,  and  pwa  is  the  reflection 
coefficient  for  radiance  coming  from  water  to  air. 

The  irradiance  just  below  the  air-water  interface,  ^(o-),  is  the  sum  of  the  direct 
sun  irradiance,  E'dn ((T ),  and  the  diffuse  irradiance  from  the  sky,  Efy  (O'): 

£,(0-)=£r(0-)+£f  (0^)  (6.8) 

Following  Diersen  and  Smith  (1997)  we  calculate  the  irradiances  just  below  the  sea 
surface  according  to  the  following  approximate  equations: 

ET{ p- )“  1-03(1  -  P,Je.  ))£?”(«" I  ( 

where  p  is  Fresnel  reflectance  of  direct  radiance  from  the  sun. 

E^(0-)~0.962Ey(0+).  (6.10) 

Here  E™(0+)  and  Ef(o+)  are  the  irradiances  just  above  the  sea  surface  of  sun  and 
sky,  respectively. 


1  -  exp  - 


K„  + 


cos  6 


\  \ 
z 

)  ) 


(6.6) 
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To  use  equation  (6.6)  we  need  to  calculate  the  apparent  optical  properties  of  water, 
Kd  and  j±d,  from  the  inherent  optical  properties  of  water.  Several  approximations  are 
available  for  these  parameters  in  the  literature.  Sathyendranath  and  Platt  (1997)  give 
the  following  equation  for  the  attenuation  coefficient  of  downwelling  irradiance: 


Kd 


a  +  b„ 


(6.11) 


where  a  is  tire  absorption  coefficient  of  water. 

This  equation  requires  knowledge  of  the  mean  cosine  of  the  downwelling  light 
stream  put.  In  the  near  sea  surface  region,  under  direct  solar  illumination,  this  parameter 
can  be  approximated  by  the  cosine  of  the  refracted  solar  zenith  angle,  Hsu,-. 


cos 

sin  1 

( i  Y 

— sin£h 

_ 

l"  JJ 

(6.12) 


where  6S  is  the  sun's  zenith  angle.  Another  extreme  approximation  for  fid  is  its 
asymptotic  value  in  sea  depth  (Haltrin,  1998): 

( 6.13) 

Z  fl 


where  fl  is  the  average  cosine  over  the  deep  water  angular  distribution  of  radiance 
(Haltrin,  1998): 


i  1/2 


a  +  3bh  +  [bb{4a  +  9bh)] 


11/2 


(6.14) 


Another  useful  approximation  for  Kd  is  given  by  Kirk  (1984): 


Kd  = 


1  +  (0.425m„ -0.19)1 


-11/2 


(6.15) 
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where  b  is  the  scattering  coefficient  of  water.  This  approximation  is  based  on 
numerical  simulations  using  the  Petzold  turbid-harbour  phase  function.  Kirk  (1991) 
also  extended  the  approximation  to  other  scattering  functions: 


Kd  = 


L 


1  +  G(jisw,g)- 


— 1 1  /  2 


(6.16) 


where 


g(p-sw>  ll. 


2.236 

8 


f  0.849  ^ 

2.447 

— 

0.739 

j 

l  8  J 

(6.17) 


Here  g  is  the  asymmetry  parameter  of  the  scattering  phase  function. 

Comparison  between  different  approximations  for  the  attenuation  coefficient  of 
downwelling  irradiance,  Kj,  is  shown  in  Figure  9.  The  inherent  optical  properties  were 
computed  using  the  one-parameter  bio-optical  model  (Section  7)  with  concentration  of 
chlorophyll  Cc=l  mg/m3  that  is  found  to  be  a  typical  value  for  Australian  coastal  water 
(Section  7).  Sun  zenith  angle  is  60  deg.  The  asymmetry  parameter  of  the  scattering 
phase  function  g=0. 95,  which  is  characteristic  value  for  seawater  (Zege,  1991). 


Figure  9.  Comparison  betxveen  different  approximations  for  the  attenuation  coefficient  of 
downwelling  irradiance.  Sun  zenith  angle  9S  =  60  deg .  Phase  function  asymmetry  parameter 
g =0.95. 
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For  the  radiance  attenuation  coefficient  we  use  the  following  expression 
(Preisendorfer,  1976): 


K-i.  ~  Dl  (a  +  bb )/ 


(6.18) 


where  D,  is  the  distribution  function.  In  the  current  work  we  assume  that  Dl  ~  1 .  This 
means  that  light  propagating  towards  the  receiver  is  lost  only  due  to  absorption  and 
backscattering.  Such  an  approximation  may  only  be  applied  for  light  reflected  from  an 
element  of  the  water  column.  For  light  reflected  from  a  target  one  should  use  the  total 
attenuation  coefficient  as  KL . 

A  more  accurate  approximation  for  Dt  is  given  by  Lee  et  al.  (1998): 


Dl  « 1.2(1  +  2.0m)0'5  ,  (6.19) 

where 

u  —  bb  /(a  +  bh ),  (  6.20) 


Equation  (6.19)  was  obtained  by  fitting  an  analytical  model  for  the  remote  sensing 
reflection  to  that  calculated  by  Hydrolight  (Mobley,  1994).  hi  further  improved 
versions  of  the  simulation  software  we  will  use  either  more  accurate  approximation 
(6.19)  or  full  numerical  solution  for  the  radiance  field  by  the  Hydrolight.  An 
investigation  in  feasibility  of  combining  the  Hydrolight  package  with  the  image 
simulation  software  is  currently  being  conducted. 

The  remote  sensing  reflection  is  defined  as  the  ratio  of  upwelling  radiance  to 
downwelling  irradiance  evaluated  just  below  the  sea  surface: 

=  £,  (0- )/£>-).  (6.21) 


From  equation  (6.6)  for  the  down-looking  sensor  (6=0)  and  deep  water  we  have: 


2  Jc(l  +  fldDL)' 


(6.22) 
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Figure  10  compares  equation  (6.22)  with  values  rff  determined  by  other 

approximations  obtained  from  numerical  calculations.  Morel  and  Gentili  (1993)  give 
the  following  expression  for  the  remote  sensing  reflectance: 
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.dp 


0.0922 

- u. 

1  —  u 


(6.23) 


Lee  et  al.  (1998)  found  the  following  approximation  from  the  Hydrolight  calculations: 


rrf  -(0.070  +  0.  155w0752)m 


(6.24) 


One  can  see  from  the  figure  that  equation  (6.22)  agrees  with  other  approximations 
only  for  low  values  of  u  which  correspond  to  open  ocean  water.  It  follows  from  one 
parameter  bio-optical  model  (Section  7)  that  the  values  of  the  parameter  u  for  most 
waters  are  less  than  0.05  (Figure  11).  Thus  the  approximation  of  (6.22)  is  mostly  useful. 


Figure  10.  Comparison  of  rff  values  as  determined  by  different  approximations. 


17 


DSTO-TR-1240 


Figure  11.  Values  of  parameter  u  -  bh  /(a  +  bb)  versus  wavelength  for  different  concentrations 
of  chlorophyll  (mg/m3). 

7.  Bio-Optical  Model 

A  bio-optical  model  relates  the  optical  properties  of  ocean  water  to  concentrations 
of  its  biological  components  such  as  chlorophyll  and  Dissolved  Organic  Matter  (DOM). 
The  model  used  in  our  simulations  is  the  one  described  by  Haltrin  (1998),  that 
summarises  empirical  data  accumulated  by  many  researchers,  mainly  for  open  ocean 
water.  For  completeness,  we  reproduce  equations  from  his  paper  in  Appendix  C. 

Here  we  compare  the  above  model  with  experimentally  measured  values  of 
absorption  and  total  attenuation  coefficient.  The  measurements  were  taken  using  a 
WetLabs  AC-9  spectral  transmissometer  at  several  different  locations  in  Australian 
coastal  waters.  First  we  consider  how  well  the  one-parameter  model  fits  the  measured 
data.  The  fitting  was  performed  by  minimising  the  sum  of  squares  of  the  relative  error 
between  the  model  and  measured  values: 


s(c>! 

/= 1 


a{X)-a, 


(7.1) 


For  the  one-parameter  model  this  sum  is  a  function  of  only  the  concentration  of 
chlorophyll  that  can  be  found  from  minimisation  procedure.  All  other  concentrations, 
i.e.  the  concentration  of  humic  and  fulvic  acids,  small  and  large  particles,  are  derived 
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from  the  concentration  of  chlorophyll  due  to  correlation  relations  between  those 
concentrations  (see  Appendix  C). 

Figure  12  shows  one  example  of  measured  data.  Here  the  absorption  (labelled  as  a 
+  wavelength)  and  total  attenuation  coefficients  (labelled  as  c  +  wavelength)  measured 
in  Sydney  Harbour  on  31  of  August  2000  are  shown.  We  can  see  a  characteristic  near¬ 
bottom  layer  of  higher  turbidity.  Having  applied  the  minimisation  procedure  at  each 
depth,  we  obtain  the  following  profiles  of  concentrations  of  various  components  of  the 
water  (Figure  13) 


Depth  (m) 


—  a412 

—  a440 
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Figure  12.  Absorption  and  total  attenuation  coefficients  versus  depth  measured  in  Sydney 
Harbour. 
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Figure  13.  Concentrations  of  water  components  derived  from  measured  absorption  and  total 
attenuation  coefficients  using  the  one-parameter  model. 


An  example  of  model  fit  to  the  measured  data  at  a  particular  depth  is  shown  in 
Figure  14.  We  can  see  that  this  particular  data  from  Sydney  Harbour  water  can  be 
described  reasonably  well  by  the  one-parameter  model.  The  maximum  error  is  18%  for 
the  absorption  value  at  555  nm.  The  fit  of  data  from  other  sites  in  Sydney  Harbour 
produced  similar  results  with  slightly  higher  maximum  error,  with  up  to  30%  for  some 
absorption  values. 
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Wavelength  (nm) 

Figure  14.  Comparison  between  the  measured  data  and  the  one-parameter  model  for  Sydney 
Harbour  water  at  depth  of  1.25  m.  The  concentration  of  chlorophyll  obtained  from  minimisation 
of  (7.1)  is  1.02  mg/m3. 

However,  the  data  from  other  Australian  coastal  waters  are  not  very  well 
described  by  the  one-parameter  model.  The  following  two  graphs  show  the 
comparison  between  the  model  and  the  data  measured  in  coastal  water  between  Cairns 
and  Thursday  Island  (Figure  15)  and  Cowley  Beach  (Figure  16).  One  can  see  that  the 
model  fit  is  not  very  good  with  errors  up  to  45%  for  absorption  values  at  555  nm.  We 
can  see  that  the  concentration  of  chlorophyll  derived  from  the  data  using  the  one- 
parameter  model  is  higher  for  these  waters.  Poor  fit  of  the  measured  data  by  the  one- 
parameter  model  suggests  that  there  is  no  strong  correlation  between  concentrations  of 
various  components  in  analysed  coastal  water. 
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Figure  15.  Comparison  between  tlie  measured  data  and  the  one-parameter  model  for  water 
between  Cairns  and  Thursday  Island  at  depth  of  15  m.  The  concentration  of  chlorophyll 
obtained  from  minimisation  of  (7.1)  is  2.65  mg/m3. 


Wavelength  (nm) 


Figure  16.  Comparison  between  the  measured  data  and  the  one-parameter  model  for  water  at 
Cowley  Beach  at  depth  of  1.05  m.  The  concentration  of  chlorophyll  obtained  from  minimisation 
of  (7.1)  is  2.89  mg/m3. 
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The  following  analysis  shows  that  in  some  cases  good  fit  still  can  be  achieved  with 
the  equations  in  Appendix  C  but  at  the  cost  of  increasing  the  number  of  independent 
parameters. 

Here  we  consider  the  data  measured  near  Kangaroo  Island,  South  Australia 
(35°39'S  and  136°00'E)  on  31  January  2000  between  17:10  and  18:15  local  time  at  depths 
of  5, 10,  and  15  metres. 

First,  we  will  try  to  find  the  concentration  of  chlorophyll  which  gives  the  best  fit  to 
the  experimental  data  using  the  least  squares  method  as  for  above  figures.  Figure  17 
shows  the  results  of  such  fitting  attempts.  The  fit  is  not  very  good  with  the  model  giving 
much  lower  values  of  scattering  coefficient,  the  difference  between  total  attenuation  and 
absorption,  than  measured. 


Wavelength  (nm) 


Figure  17.  Comparison  betiveen  the  measured  data  and  the  one-parameter  model  for  water  off 
Kangaroo  Island  at  depth  of  5  m.  The  concentration  of  chlorophyll  obtained  from  minimisation 
of  (7.1)  is  0.212  mg/m3. 

To  achieve  a  better  fit  of  both  coefficients  we  need  to  use  more  than  one 
independent  parameter.  For  this  purpose  we  suggest  using  experimental  values  of 
absorption  coefficient  and  equations  (C.l),  (C.8),  and  (C.9)  to  find  the  concentration  of 
chlorophyll  which  gives  the  best  fit  of  experimental  data  using  equation  (C.l).  Then  we 
assume  the  concentrations  of  small  Cs  and  large  C/  particles  to  be  independent  and  find 
them  by  fitting  experimental  values  of  scattering  coefficient  by  equation  (C.2)  using  the 
least  squares  method.  Using  experimental  values  from  17:10  and  depth  of  5  m  we  find 
the  following  values  of  concentrations:  Cc=0.123  mg/m3,  0=0.038  g/m3,  and 
0=0.243  g/m3.  Figure  18  shows  the  results  of  this  fitting  procedure. 
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Wavelength  (nm) 

Figure  18.  Three-parameter  model  fitting.  Cc=0.123  mg/m3,  C s=0.038  g/m3,  and  C\=0.243 
g/m3. 


Figure  19.  Backscattering  coefficient.  One-parameter  model:  Cc=0.212  mg/m3,  Cs=0.0038 
g/m3,  and  Q =0.163  g/m3.  Three-parameter  model:  Cc=0.123  mg/m3,  Cs=0.038  g/m3,  and 
Q =0.243  g/ m3. 

Figure  19  shows  that  the  backscattering  coefficient  for  the  three-parameter  model 
is  higher  than  for  the  one-parameter  model.  This  is  because  the  relative  concentration 
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of  small  particles  is  higher  for  the  three-parameter  model  than  for  the  one-parameter 
model.  The  Kopelevich  hydrosol  parameter  determines  the  relative  volume  of  small 
particles  and  can  be  calculated  as  Cs  /(Cs  +2C,).  For  the  one-parameter  model  and  for 

the  data  analysed,  this  parameter  is  0.01  while  for  the  three-parameter  model  the  value 
of  Kopelevich  hydrosol  parameter  is  0.07,  clearly  indicating  a  higher  relative 
concentration  of  the  small  particles.  The  asymmetry  parameter  of  the  Kopelevich  phase 
function  for  the  concentration  of  particles  obtained  from  the  one-parameter  model  is 
g=0.953.  Concentrations  of  particles  obtained  from  the  three-parameter  model  give  a 
value  of  asymmetry  parameter  of  g=0.934.  The  Kopelevich  phase  function  derived 
from  the  transmissometer  data  for  water  off  Kangaroo  Island  using  one  and  three- 
parameter  model  is  presented  in  Figure  20. 


Figure  20.  Kopelevich  phase  function  derived  from  tire  AC-9  data  for  xvater  off  Kangaroo  Island. 
One-parameter  model:  Cc-0.212  mg/m3,  Cs=0.0038  g/m3,  and  Ci =0.163  g/m3,  g=0.953. 
Three-parameter  model:  Cc-0.123  mg/m3,  C s=0.038  g/ m3,  and  Q=0.243  g/ m3,  g=0.934. 


8.  Target  Rendering 

Target  rendering  is  based  on  the  ray  tracing  technique.  For  each  ray  we  determine 
whether  it  crosses  the  target  surface.  If  it  does,  we  find  the  coordinates  of  the 
intersection,  illumination  of  the  surface  element,  and  reflected  radiance  towards  the 
receiver. 
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A  ray  can  be  described  by  the  following  set  of  equations: 
x  =  x0  +  vxt, 

h  =  ho+V’ 

(8.1) 

z  =  z0+  vj,  v 

v2  +  V2  +  V2  =  1. 

x  y  z 

For  the  ray  in  air  (x0  ,y0,z0)  are  the  coordinates  of  the  focal  point  of  the  camera.  For 
the  underwater  part  of  the  ray  (x0,y0,z0)  are  the  coordinates  of  the  ray  intersection 
with  the  sea  surface.  Here  v,,,v  ,andv_  are  the  directional  cosines  of  the  ray,  t  is  the 
distance  from  the  ray  starting  point  to  the  current  point  (x,  y,  z). 

The  target  surface  can  be  described  as 

F(x,y,z)=  0  (8.2) 


Substituting  equation  (8.1)  into  (8.2)  we  have  an  equation  for  the  ray  parameter  t .  For 
simple  geometrical  shapes  such  as  spheres,  cylinders,  ellipsoids  etc.  this  equation  can 
be  solved  analytically.  If  the  equation  has  multiple  roots,  we  take  the  one  with  the 
minimum  value  of  the  parameter  t .  We  can  then  find  the  unit  vector  n  normal  to  the 
target  surface  at  the  intersection  point  and  calculate  irradiance  at  the  surface  element: 

=  J|(n-sJZ(A,s)fK2,  (8.3) 

2  K 


where  L( s)  is  the  underwater  radiance  in  the  direction  specified  by  the  unit  vector  s, 
dQ.  is  the  differential  solid  angle.  The  integration  in  equation  (8.3)  is  performed  over 
the  hemisphere  outside  the  target. 

We  currently  assume  that  the  target  surface  is  lambertian.  Therefore,  the  radiance 
reflected  towards  the  receiver  will  be 

Lr{X)  =  ^-R(X),  (8.4) 

Jl 

where  f?(A)  is  the  spectral  reflectance  of  the  target  surface. 

A  complex  three-dimensional  target  can  be  formed  by  combining  several  simple 
elements.  In  this  case  we  need  to  solve  the  system  of  equations  (8.1)  and  (8.2)  for  each 


26 


DSTO-TR-1240 


element  of  the  target  and  take  the  solution  with  the  minimum  value  of  the  ray 
parameter  t .  The  sea  bottom  can  be  thought  of  as  another  element  of  the  target.  In 
most  cases  it  will  be  an  infinitely  wide  plane. 

Target  blur  caused  by  light  scattering  in  seawater  is  currently  modelled  using  a 
simple  empirical  Point  Spread  Function  (PSF)  of  ocean  water  (Voss,  1991).  The  PSF  is 
applied  in  postprocessing  to  the  'underwater'  part  of  the  simulated  image,  i.e.  to  the 
radiance  coming  from  below  the  sea  surface.  To  do  this  we  save  separately  the  radiance 
from  underwater  L(™]  and  the  radiance  reflected  from  the  sea  surface  ,  where  m 

and  n  are  the  pixel  indices,  and  then  convolve  L(™]  with  the  PSF: 


ni'=~r  n'=~r 


(8.5) 


where  Tmn  is  the  underwater  path  of  the  ray  from  pixel  (m,n),  is  the  angle 

between  underwater  paths  of  the  rays  from  pixels  (m,n)  and  (m',n'): 


e"=ocos(v 


mn  m-m' ,n-ri  .  mn  m~m' ,n~ri  .  mn  m-m 
'  V  y  '  -  - 


(8.6) 


where  (v™ ,  v"w ,  v"" )  are  the  directional  cosines  of  the  underwater  path  of  the  ray  from 
pixel  (ra,n). 

We  can  see  from  equations  (8.5)  and  (8.6)  that  to  perform  such  a  convolution  we 
need  to  also  save  the  underwater  paths  and  directional  cosines  for  each  ray.  The 
convolution  itself  is  rather  a  slow  process.  The  use  of  the  Modulation  Transfer  Function 
(MTF)  and  the  Fast  Fourier  Transform  (FFT)  would  increase  the  speed  of  this  process 
significantly.  Flowever,  the  MTF  is  not  suitable,  for  example,  in  the  case  of  a  three- 
dimensional  target  with  a  significant  difference  in  range  from  the  camera  to  different 
parts  of  the  target.  This  is  because  the  MTF  technique  does  not  allow  variable  blur 
across  the  image.  In  this  case  the  use  of  PSF  is  preferable  because  PSF  explicitly 
includes  the  range  to  the  various  points  of  the  target  through  the  ray  optical  path  T  . 

9.  Sensor  model 


All  of  the  above  steps  for  calculating  the  radiance  reflected  or  backscattered  from 
various  parts  of  the  simulated  scene  eventually  produce  the  spectral  radiance  arriving 
at  the  receiver.  To  simulate  the  image  generated  by  a  specific  sensor  we  need  to 
convolve  the  radiance  at  the  receiver  with  the  spectral  sensitivity  of  the  sensor. 
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In  the  case  of  the  human  eye  or  digital  camera  the  spectrum  must  be  transformed 
into  three  RGB  values  to  produce  an  image: 


4(A)® 
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where  Lr  is  the  radiance  at  the  receiver,  x,y,z  are  the  colour-matching  functions, 
X,Y,Z  are  the  three  standard  primaries  defined  by  International  Commission  of 
Illumination  (CIE)  in  1931,  and  the  operator  ®  denotes  convolution.  The  value  Y  also 
represents  the  luminance  that  will  be  used  in  the  next  section  for  comparison  with  the 
measurements. 

In  the  case  of  a  hyperspectral  sensor  there  are  several  wavebands,  each  having  a 
specific  bandwidth.  Average  radiances  within  these  wavebands  can  be  calculated  as 
follows: 
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where  xt  (A)  is  the  spectral  sensitivity  function  of  bands. 


10.  Simulation  Examples 


Now  we  consider  several  examples  of  simulated  images  for  different  scenarios. 
Figure  21  shows  a  comparison  of  a  simulated  image  (left)  with  a  real  photograph  (right). 
The  photograph  was  taken  with  a  Nikon  D1  digital  camera  in  the  Gulf  of  St  Vincent  on 
20th  of  December  2000  at  1 1  a.m.  Field  of  view  of  the  camera  lens  was  51  degrees.  The 
sun  elevation  above  horizon  was  59  degrees.  The  camera  was  pointed  directly  away 
from  the  sun.  Estimated  wind  speed  was  about  5  m/s.  The  image  on  the  left  in  Figure  21 
is  simulated  for  the  same  sun  position  with  respect  to  the  camera  and  the  same  field  of 
view.  Sea  surface  realisation  is  calculated  for  wind  speed  of  5  m/s  using  Apel's  wave 
spectrum  (Apel,  1994).  Sky  radiance  is  calculated  from  the  model  of  Zibordi  and 
Voss  (1989). 

The  optical  properties  of  seawater  are  calculated  from  the  one-parameter  bio-optical 
model  of  Haltrin  (1998)  with  a  concentration  of  chlorophyll  value  of  1.4  mg/m3.  This 
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value  is  inferred  from  fitting  the  one-parameter  model  to  the  data  measured  on  the  day 
with  the  AC-9  spectral  transmissometer. 


Figure  21.  Comparison  of  simulated  image  (left)  with  real  photograph  (right). 


Qualitatively,  the  simulated  image  compares  reasonably  well  with  the  photograph 
especially  in  the  sky  area.  The  sea  surface  looks  somewhat  darker  in  the  simulated 
image.  One  of  the  reasons  for  this  may  be  as  follows.  The  model  always  assumes  the 
presence  of  small  capillary  waves  in  every  pixel  footprint  at  the  sea  surface.  In  reality, 
however,  this  may  not  be  the  case  and  some  waves  may  be  rather  smooth  on  the  scale 
of  an  individual  pixel  footprint.  This  may  result  in  higher  reflection  of  the  sky  radiance 
from  the  sea  surface. 

With  respect  to  quantitative  comparison.  Figure  22  compares  luminance  calculated 
during  the  image  simulation  with  the  measurements  made  at  the  time  of  taking  the 
photograph.  The  solid  curve  represents  the  luminance  calculated  for  the  vertical  centre¬ 
line  of  the  image  with  some  extension  beyond  the  upper  boundary  of  the  image.  Zero 
angle  in  Figure  22  corresponds  to  the  horizon  with  the  negative  angles  representing 
atmosphere  and  the  positive  angles  representing  the  sea  surface.  The  measurements 
were  made  with  a  Tektronix  LumaColor  photometer  J17  using  luminance  sensor-head 
J1803  with  a  field  of  view  of  8  degrees.  The  sensor  head  was  hand  held  during  the 
measurements,  so  the  look  angle  of  the  sensor  and,  therefore,  position  of  the  bars  in 
Figure  22  is  rather  approximate.  Taking  this  uncertainty  into  account,  the  calculated 
luminance  compares  reasonably  well  with  the  measured  data. 
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Figure  22.  Comparison  of  simulated  luminance  (solid  curve)  with  measurements  (red  bars).  The 
width  of  red  bars  represents  a  field  of  view  of  8  degrees. 


Figure  23  demonstrates  the  model's  capability  to  render  three-dimensional 
underwater  objects.  In  the  simulated  image  we  see  a  red  cylindrical  object,  three  grey 
spheres  of  various  sizes,  and  a  model  of  a  Manta  mine.  The  cylindrical  object  is 
composed  of  three  basic  geometrical  objects:  a  cylinder  and  a  disc  at  each  end  of  the 
cylinder.  The  model  of  the  Manta  mine  is  composed  of  a  red  truncated  cone  with  a  red 
disc  on  top  of  the  cone  and  a  black  disc  on  top  of  the  red  disc.  The  water  depth  is  2.5  m. 
The  optical  properties  of  water  are  calculated  from  the  one-parameter  model  of  Haltrin 
(1998)  with  a  concentration  of  chlorophyll  value  of  1  mg/ m3.  The  bottom  is  sand  with 
reflectivity  changing  linearly  from  40%  at  400  nm  to  55%  at  700  nm.  The  image  is 
disturbed  by  a  wavy  sea  surface  that  is  generated  from  Apel's  wave  spectrum  using  a 
wind  speed  value  of  0.5  m/ s.  The  wind  direction  is  at  45  deg  to  the  look  direction.  The 
sun  elevation  is  26  deg  and  its  azimuth  is  20  deg  to  the  right  from  the  look  direction. 
The  sensor  height  above  sea  surface  is  5  m,  zenith  look  angle  is  30  degrees,  and  sensor 
field  of  view  is  24  degrees.  The  blur  due  to  scattering  is  modelled  by  applying  the 
point-spread  function  to  the  underwater  part  of  the  image  as  described  above  in 
Section  8. 
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Sun 


Figure  23.  Simulated  image  ofundenoater  objects. 

Another  passive  system  that  can  be  simulated  with  our  model  is  a  hyperspectral 
sensor.  Airborne  hyperspectral  sensors  are  predominantly  "pushbroom"  type  sensors. 
In  one  swath  the  sensor  takes  an  image  of  a  narrow  strip  of  water  across  the  direction 
of  aircraft  flight.  The  image  is  built  up  from  the  cross-track  lines  of  pixels  while  the 
aircraft  moves  along.  The  image  is  taken  in  several  narrow  spectral  bands,  the 
appropriate  choice  and  combination  of  which  may  allow  for  the  detection  of 
underwater  objects. 

The  input  to  the  model  in  this  case  consists  of  parameters  such  as  across-track 
field-of-view,  number  of  across-track  pixels,  aircraft  altitude,  speed  and  headings,  and 
spectral  bands.  Environmental  parameters  include  wind  speed,  parameters  of  the  bio- 
optical  model  or  the  spectral  optical  properties  of  water,  sun  altitude  and  azimuth, 
bottom  depth  and  profile.  The  other  parameters  include  target  shape,  reflectance,  and 
position  with  respect  to  image  centre.  One  of  the  output  files  contains  spectral 
radiances  at  specified  spectral  bands,  interleaved  by  pixel.  The  file  can  be  read  by 
standard  processing  software  such  as  ENVI. 

Figure  24  shows  an  example  of  a  simulated  image  and  the  result  of  simple 
processing  in  ENVI  using  the  ratio  of  the  band  at  530  nm  (green)  to  the  band  at  470  nm 
(blue).  The  direction  of  flight  is  from  the  right-hand  side  of  the  image  to  the  left,  wind 
is  5m/s  at  45  degrees  to  the  direction  of  flight,  with  the  sun  at  25  degrees  above 
horizon  right  behind  the  aircraft.  The  bottom  is  sand  with  reflectivity  varying  linearly 
from  40%  at  400  nm  to  55%  at  700  nm.  Bottom  depth  varies  from  8  m  at  the  right  end  to 
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9.5  m  at  the  left  end.  The  aircraft  altitude  is  500  m,  cross-track  field  of  view  is  40 
degrees,  resulting  in  a  364  m  swath.  There  is  a  target  in  the  middle  of  the  surveyed 
area.  The  target  is  a  grey  rectangle  25m  x  25m  with  45%  reflectance. 


b) 

Figure  24.  Example  of  simulated  hyperspectral  image,  a)  RGB  "true  colour "  image;  b)  processed 
in  ENVI  using  band  ratio. 


In  this  simulated  case  even  the  simple  analysis  allows  us  to  see  a  target  which  is 
not  visible  in  the  conventional  image.  The  processing  also  demonstrates  that  the 
shallower  water  the  greener  it  appears  which  is  in  agreement  with  everyday 
experience. 


11.  Conclusion 


A  computer  model  has  been  developed  to  simulate  the  performance  of  passive 
optical  imaging  sensors,  such  as  a  human  eye,  a  conventional  digital  camera  or  a 
hyperspectral  sensor.  The  model  is  based  on  the  ray  tracing  technique  and  implements 
several  basic  models  describing  various  elements  of  image  formation.  These  models  are 
the  analytical  model  of  sun  and  sky  radiance  by  Zibordi  and  Voss  (1989)  and  the  model 
of  wind  driven  waves  on  sea  surface  based  on  the  directional  wave  spectrum  by 
Apel  (1994).  The  sun  and  sky  reflection  by  capillary  waves  on  the  sea  surface  is  based 
on  the  model  of  Zeisse  (1995).  Optical  properties  of  seawater  are  modelled  by  the  one- 
parameter  bio-optical  model  of  Haltrin  (1998).  Several  empirical  and  semi-empirical 
relations  are  used  for  describing  such  elements  as  radiance  and  irradiance  transfer 
through  a  wavy  sea  surface,  and  image  blur  due  to  light  scattering  in  ocean  water.  An 
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analytical  expression  for  radiance  scattered  into  a  ray  path  has  been  developed  in  the 
present  work. 

All  these  models  are  based  on  approximations  of  various  degrees  of  accuracy. 
Some  of  them  were  verified  in  the  present  work  against  other  models  and  ground  truth 
measurements.  The  model  produces  realistic  images  of  sea  surface  and  underwater 
targets  but  further  verification  of  the  simulation  model  and  its  components  is  required 
and  is  planned  in  the  nearest  future. 

The  future  work  will  include  comparison  of  simulations  with  hyperspectral  data 
and  ground  truth  measurements.  Further  improvements  to  the  model  elements  will  be 
made.  It  is  planned  to  expand  the  model  to  incorporate  simulation  of  the  performance 
of  active  imaging  systems  such  as  Laser  Range  Gated  (LRG)  cameras  and  the  Streak 
Tube  Imaging  Lidar  (STIL).  Also,  atmospheric  effects  on  light  propagation  near  the  sea 
surface  will  be  included  in  the  model. 
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Appendix  A:  Analytical  Model  of  Sky  Radiance 

The  following  model  for  sky  radiance  is  essentially  from  Zibordi  and  Voss  (1989). 
Solar  radiance  just  above  the  sea  surface  can  be  represented  as  the  sum  of  two 
components:  direct  and  diffuse  radiance.  The  direct  solar  irradiance  can  be  calculated 
using 

E(X)  =  E0D(X)exp(~T(X)m),  (A.l) 

where  E0D(X)  is  the  extraterrestrial  solar  irradiance  corrected  by  the  Sun-Earth 
distance: 


E,D(l)=E„(X 


1  +  ecos 


2ft(D-3)V 
365  . 


(A.2) 


with  E0  (A)  being  the  extraterrestrial  solar  irradiance  at  wavelength  Z,£  eccentricity  of 
the  Earth  orbit,  and  D  the  Julian  day.  In  equation  (A.l)  t(A)  is  the  total  optical 
thickness  of  the  atmosphere,  and  m  is  the  relative  air  mass  (Kasten,  1966),  which  may 
be  approximated  as  l/cos0o  for  solar  zenith  angles  (60)  less  than  60°.  For  unit  air 
mass. 


t(Z)—  Ty?(A)+T/i(A)+T0j  (X)+tw(X)+tg(X) 


(A.3) 


Various  optical  thicknesses  in  the  right  hand  side  of  equation  (A.3)  can  be 
computed  by  various  empirical  equations  as  follows. 

The  scattering  optical  thickness  of  permanent  gases  (molecular,  or  Rayleigh, 
optical  thickness)  can  be  computed  (Marggraf  and  Griggs,  1969)  as 

tR  (^) =  0.0088A-4 15+0  2A ,  (A.4) 


where  X  is  in  units  of  micrometers. 

The  absorption  and  scattering  optical  thickness  of  aerosol,  assuming  a  radius 
power  law  distribution  of  particles,  can  be  given  by  (Angstrom,  1961) 

IA(X)  =  aX~pr  (A.5) 

where  a  and  respectively  relate  to  the  concentration  and  size  distribution  of 
particles.  Typical  values  of  the  wavelength  coefficient  ft  lie  in  the  range  0.8  -  2. 
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Turbidity  coefficient  a  can  be  calculated  as  a  =  a00.5p  13 .  Typical  values  of  the 
Angstrom  turbidity  coefficient  (X0  are  given  in  Table  1. 


Table  1.  Typical  Angstrom  turbidity  coefficients  (Leckner,  1978) 


Latitude 

60°N 

45°N 

30°N 

00 

Low  limit 

0.01 

0.047 

0.047 

0.047 

Medium 

0.057 

0.077 

0.097 

0.117 

High  limit 

0.093 

0.187 

0.375 

0.375 

The  absorption  optical  thickness  of  ozone  can  be  computed  through 

*o3W=Q>3£03W,  (A.6) 

where  K0}  (A)  (cm-1)  is  the  spectral  absorption  coefficient  of  the  ozone  (Table  2, 
Vigroux,  1953)  and  C0  (cm)  is  the  ozone  concentration. 


Table  2.  Spectral  absoiption  coefficient  of  ozone  (Vigroux, 1953). 


A 

—ffo. 

A 

Kp. 

A 

K0. 

A 

Ko- 

A 

Kn 

290 

38 

445 

0.003 

515 

0.045 

585 

0.118 

700 

0.023 

295 

20 

450 

0.003 

520 

0.048 

590 

0.115 

710 

0.018 

300 

10 

455 

0.004 

525 

0.057 

595 

0.12 

720 

0.014 

305 

4.8 

460 

0.006 

530 

0.063 

600 

0.125 

730 

0.011 

310 

2.7 

465 

0.008 

535 

0.07 

605 

0.13 

740 

0.01 

315 

1.35 

470 

0.009 

540 

0.075 

610 

0.12 

750 

0.009 

320 

0.8 

475 

0.012 

545 

0.08 

620 

0.105 

760 

0.007 

325 

0.38 

480 

0.014 

550 

0.085 

630 

0.09 

770 

0.004 

330 

0.16 

485 

0.017 

555 

0.095 

640 

0.079 

780 

0 

335 

0.075 

490 

0.021 

560 

0.103 

650 

0.067 

790 

0 

340 

0.04 

495 

0.025 

565 

0.11 

660 

0.057 

800 

0 

345 

0.019 

500 

0.03 

570 

0.12 

670 

0.048 

810 

0 

350 

0.007 

505 

0.035 

575 

0.122 

680 

0.036 

820 

0 

355 

0 

510 

0.04 

580 

0.12 

690 

0.028 

830 

0 

The  absorption  optical  thickness  of  water  vapour  (Iqbal,  1983)  can  be  calculated  by 


0.2385 Kw(k)Cw 
(1  +  20.07 Kw{X)CwyA5 


(A -7) 
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where  Kw  (A)  (cm-1)  is  the  spectral  absorption  coefficient  of  water  vapour  (Leckner, 
1978)  and  Cw  (cm)  is  the  precipitable  water  of  the  atmosphere. 

The  absorption  optical  thickness  of  the  so-called  uniformly-mixed  gases  (mainly 
oxygen  and  carbon  dioxide)  can  be  calculated  (Iqbal,  1983)  by 


T  iaikg&) 

gK  '  [1  +  118.93A:g(A)]045 


(A.8) 


where  KG  (A)  is  the  effective  absorption  coefficient  (Leckner,  1978).  In  the  visual  part 
of  the  spectrum  the  water  vapour  absorption  and  the  uniformly-mixed  gases 

absorption  can  be  neglected. 

The  diffuse  part  of  the  solar  radiance  is  due  to  scattering  of  sun  light  by  particles 
and  molecules  in  the  atmosphere.  For  computation  of  the  diffuse  radiance  it  is 
necessary  to  separate  the  total  aerosol  optical  thickness  T  A  (A)  into  its  scattering  T^(A) 

and  absorption  tAa  (A)  components: 

?<(A)  =  T^(A)+tJA).  (A.9) 

This  can  be  done  through  the  aerosol  single  scattering  albedo  C0A  (A),  which  is  the  ratio 
of  energy  scattered  to  total  attenuation  by  aerosol,  and  defined  by 

T A  vU 

so  that 

tA:(a)  =  o>,(A)tJ/.),  (A.ll) 

T^W=[1  (A.12) 

(0A  (A)  strongly  depends  on  the  optical  properties  of  the  aerosol  particles.  Generally  it 
is  assumed  invariant  with  wavelength  in  the  visible  and  near  infrared  and  ranges 
between  0.6  and  1.0,  with  (0A  ~  0.6  for  urban/ industrial  aerosol,  0)A  ~  0.9  for 
continental  aerosol  and  0)A  =1.0  (pure  scattering)  for  maritime  aerosol  (Zege  et  al, 
1991). 

The  phase  function  defines  the  angular  distribution  of  sky  radiance  due  to 
scattering  of  diffusing  constituents.  The  molecular  component  (Rayleigh  scattering)  is 
given  by  the  classical  relationship  for  isotropic  scattering 
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/>#(4')=|(l  +  cos2‘P),  (A.13) 

where  *F  is  the  scattering  angle. 

The  aerosol  phase  function  can  be  approximated  by  the  Two  Term  Henyey- 
Greenstein  (TTHG)  analytical  function 

r+7  (a.m) 

(1  +  gi  ~2g,  cos^)  (l  +  g;+2g2cosxFj 

where  parameters  a,  gi,  and  g2  (asymmetry  factors)  depend  on  the  aerosol  size 
distribution  and  wavelength,  however,  the  latter  is  less  significant.  The  asymmetry 
parameters  for  the  TTHG  function  (a=0.962,  gi=0.713,  and  g2=0.759)  were  taken  from 
Kattawar  (1975).  These  values  were  successfully  used  in  remote  sensing  applications 
with  clear  sky  conditions  (Gordon  et  ah,  1983). 

The  global  phase  function  of  the  atmosphere  is  given  by 
„  W  =  T,(A)f„('P)+r„(A)P,(y) 

T„  (/.)+  ZA.  ()■-) 

with  the  normalisation  condition 
1  K 

— J  F(vF)sin'F<AF  =  1, 

2  o 

and  the  scattering  angle  *F  for  sky  radiance  defined  by 
'F  =arccos(cos0cos0o  -  sin  0  sin  60  cos((f)  -  0O )) .  (A.17) 

where  60,(j)0  and  9,(j)  are  the  sun  and  sensor  zenith  and  azimuth,  respectively. 

Calculation  of  sky  radiance  is  based  on  the  approximate  solution  of  the  radiative 
transfer  equation  given  by  Sobolev  (1975)  for  a  plane  parallel,  homogeneous  and 
cloudless  atmosphere  bounded  by  a  lambertian  surface.  This  solution  allows  the  exact 
computation  of  the  first  order  of  scattering,  while  multiple  orders  are  accounted  for 
through  the  phase  function  averaged  over  the  azimuth.  Absorbing  constituents  are 
assumed  condensed  in  a  thick  layer  at  the  top  of  the  atmosphere,  which  enables 
independent  computation  of  absorbing  and  scattering  processes.  In  such  a  situation  the 
sky  radiance  can  be  given  by 


(A.15) 


(A.16) 
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Lsky  (A)  =  cr(rrf  (A))cos0o  exp1 


K 


-Tgfa) 

COSfL 


(A.18) 


where  rrf(A)  and  tb(A)  are  the  total  optical  thickness  for  scattering  and  absorbing 
constituents,  respectively,  and  are  given  by 


?rf(A)  =  T*(A)+T^(A), 

tu  (A)  =  t0  (A)+th,  (A)+  tg  (A)+t,0  (A) . 


(A.19) 

(A.20) 


The  term  <j(?d  (A))  is  the  scattering  transmission  coefficient  of  the  atmosphere  and, 
following  Sobolev  (1975),  is  given  by 

/  / *  w  fll  -  p(AMA,0)+  2p(A)}f?(A,e0) 

=  4  +  (3-P1)Ll-p(A)K(A) 


f 


exp 


~Td  (a) 


+  exp- 


■  rd  (A) 


v 


j 


cos  6  x  cos  90 

+[pM-  (3  +  P\  )cos0  coseioiT,  (t,  (Aft 


(A.21) 


where  the  term  ox  (rd  (A))  expresses  the  single  scattering  transmission 
■T,(A)  -T^A) 


<Ji(rrf(A))= 


exp 

1  COS0 


-  exp- 


cos  0 


COS0  -COS0n 


(0*0O). 


(A.22) 


Px  is  the  first  term  in  the  Legendre  polynomial  expansion  of  the  phase  function 


Px  =  —  |  P(xP)cos  'F  sin  WcN? , 


R( A,  £  )  (with  £  =  9, d0  )  is  an  auxiliary  function 


(A.23) 


R( A,£)  =  1  +  -cos£  + 


1-— cos  £ 


iexp| 


-*A A) 


COS^ 


(A.24) 
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and  p(/ 1)  is  spectral  albedo  of  the  surface  (assumed  lambertian). 
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Appendix  B:  Wavenumber  Spectrum 


Here  we  present  a  short  summary  of  the  Apel's  directional  wavenumber  spectrum 
following  Apel  (1994).  The  spectrum  can  be  written  as  follows: 

¥(*)=¥ v{k)D{k,(f>)r  (B.l) 

where  k  is  the  wavenumber  and  0  is  the  azimuth  angle.  y/(k)  is  the  portion  of  the 
spectrum  that  is  independent  on  the  azimuth  angle,  0.  It  can  be  represented  as  a 
multiplication  of  the  five  factors: 

¥(k)=A-L0-Jp-k^-H,  (B.2) 

Here,  the  spectral  constant 

A  =  0.00195,  (B.3) 

the  low-frequency  roll-off 

L0  =exp[-(kp/k}\,  (B.4) 

where 


g 

V2 u2 


(B.5) 


with  g  being  the  specific  gravity  and  u  the  wind  speed  measured  at  ten  metres  height 
above  sea  surface. 

The  spectral  peak 


(B.6) 


For  fully  developed  seas,  the  base,  T,  is  1.7,  and  the  width  parameter  of  the  spectral 
peak,  8,  is  0.40. 

A  combination  of  high  frequency  effects  is  given  by 

(R7> 
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where 

Rm  = - 7— - -r-,  with  krn  =100  rad/m. 

\  +  (k/kJ 

Rre,  =  aksech[(k  -  kres )/  kM],  with  a= 0.8;  lcrcs=400;  kw= 450. 

S(u)=  expfs,  +  s2  (l  -  e~u'u"  )Jlne  1  o}, 

where  5,  =  -4.95,  s2  =3.45,  and u  =  4.7 m/s 

K,  =  exp(-  k2 1 kL  \  kd,  *  6283. 

The  directional  part  of  the  spectrum  is: 

£>(M)  =  ex p\~4frr-  ■ 

L  20;  \k) 

2<P;(k)=[i.U  +  5.o(klkpY2]' , 

where  <pu  is  the  wind  direction. 


(B.8) 

(B.9) 

(B.10) 

(B.ll) 

(B.12) 

(B.13) 


44 


DSTO-TR-1240 


Appendix  C:  Bio-Optical  Model 


This  Appendix  presents  a  summary  of  the  one-parameter  bio-optical  model  of 
Haltrin  (1999).  Absorption  is  described  as  the  sum  of  absorption  coefficients  of  pure 
water,  chlorophyll,  and  two  components  of  Dissolved  Organic  Matter  (DOM):  fulvic 
and  humic  acids: 

a{X)  =  aw  (A)+  0.06ac°  (A)(cc  / Cc° )° 6°"  +  a°fCf  exp(-  kfX)+  a°hCh  exp(-  khX), 

(C.1) 

where  aw(X)  is  the  absorption  coefficient  of  pure  water  in  mA  (Pope  and  Fry,  1997),  X 
is  the  vacuum  wavelength  in  nm,  a°c(X)  is  the  specific  absorption  coefficient  of 
chlorophyll  in  niA  (Prieur  and  Sathyendranath,  1981),  Cc  is  the  total  concentration  of 
chlorophyll  in  mg/m3  (C°  =1  mg/m3),  a°f  =  35.959m2  / mg  is  the  specific  absorption 
coefficient  of  fulvic  acid,  kf  =  0.0189«77i_1 ;  a°h  =  18.828m2 /mg  is  the  specific 
absorption  coefficient  of  humic  acid,  kh  =0.01105 nm~l;  Cf  and  Ch  are,  respectively, 
concentrations  of  fulvic  and  humic  acids  in  mg/m 3  (Carder  et  al,  1989). 

The  scattering  coefficient  b(X)  and  backscattering  coefficient  bB  (/l)  are  calculated 
according  to  following  equations  (Haltrin,  1999): 

4(q  =  »„W+*,°WC,  +i,"U)C„  (C.2) 

(A)  =  0.56.  (A)+  B,b°  (A)C,  +  B,bf  (A)C, ,  (C.3) 

here 


It  K 

Bs  =  0.5  J ps(d)sin6d9  =  0.039,  B,  =  0.5  jp,{d)sm6de  =  6.4-1 0~4 ,  (C.4) 

n! 2  k!2 


Bs  is  the  probability  of  backscattering  by  small  particles,  Bi  is  the  probability  of 
backscattering  by  large  particles,  bw(X)  is  the  scattering  coefficient  by  pure  water  in 

m~x ,  65°(/l)  and  tf{X)  are,  respectively,  the  specific  scattering  coefficients  in  m2/g  for 
small  and  large  particulate  matter,  Cs  and  C/  are,  respectively,  concentrations  in  g/m3  of 
small  and  large  particles.  The  following  equations  are  given  by  Haltrin  (1999): 


bw{h)=  (5.826 -lO^m-1 


400 


\4.322 

) 


(C.5) 
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6°(/l)=(l.l513m2/*)| 


(  400" 
V  ^  > 


1.7 


b°  (A)  =  (o.341  lm2  / g)\ 


400 


\0.3 


(C.6) 


(C.7) 


Haltrin  (1999)  finds,  from  experimental  data  of  other  researchers,  that  all  five 
concentrations  in  above  equations  correlate  with  each  other  and  suggests  using  the 
concentration  of  chlorophyll  as  the  main  parameter  and  express  other  concentrations  as 


follows: 

Cf  =  1.74098  •  Cc  •  exp(0. 12327  •  Cc )  (C.8) 

Ch  =  0.19334  •  Cc  •  exp(0. 12343  •  Cc )  (C.9) 

Cv  =  0.01739  -Cc-exp(0.1 163 1-Cc)  (C.10) 

C,  =  0.76284  ■  Cc  •  exp(0. 03092  •  Cc )  (C.ll) 


46 


DISTRIBUTION  LIST 


Simulation  of  Passive  Imaging  of  Underwater  Objects 
Alexei  Kouzoubov 

AUSTRALIA 


DEFENCE  ORGANISATION 

Task  Sponsor  DNW  2  copies 

S&T  Program 

Chief  Defence  Scientist  ] 

FAS  Science  Policy  \  shared  copy 

AS  Science  Corporate  Management 
Director  General  Science  Policy  Development 
Counsellor  Defence  Science,  London  (Doc  Data  Sheet) 

Counsellor  Defence  Science,  Washington  (Doc  Data  Sheet) 

Scientific  Adviser  Joint 
Navy  Scientific  Adviser 

Scientific  Adviser  -  Army  (Doc  Data  Sheet  and  distribution  list  only) 

Air  Force  Scientific  Adviser  (Doc  Data  Sheet  and  distribution  list  only) 

Director  Trials  (Doc  Data  Sheet  and  distribution  list  only) 

Aeronautical  and  Maritime  Research  Laboratory 
Director 

Chief  of  Maritime  Operations  Division 
Research  Leader,  Mine  Warfare 
Research  Leader,  Maritime  Sensor  Systems 
Head  Maritime  Optronics 

Author(s):  A.Kouzoubov  2  copies 

V.K.Shettigara 
S.Sutherland 
P.  Chappie 
P.Mulheam 

DSTO  Library 
Library  Fishermans  Bend 
Library  Maribymong 
Library  Edinburgh 
Australian  Archives 
Library,  MOD,  Pyrmont 
Library,  MOD,  HMAS  Stirling 
*US  Defense  Technical  Information  Centre, 

*UK  Defence  Research  Information  Center, 

^Canada  Defence  Scientific  Information  Service, 


(Doc  Data  sheet  only) 
(Doc  Data  sheet  only) 


2 

2 

1 


‘'NZ  Defence  Information  Centre, 
National  Library  of  Australia 


1 


Capability  Development  Division 

Director  General  Aerospace  Development  (Doc  Data  Sheet  only) 
Director  General  Maritime  Development 


Knowledge  Staff 

Director  General  Command,  Control,  Communications  and  Computers  (DGC4) 
(Doc  Data  Sheet  only) 


Navy 

SO  (Science),  Director  of  Naval  Warfare,  Maritime  Headquarters  Annex, 
Garden  Island,  NSW  2000.  (Doc  Data  Sheet  only) 

Intelligence  Program 

DGSTA  Defence  Intelligence  Organisation 

Manager,  Information  Centre,  Defence  Intelligence  Organisation 

Corporate  Support  Program 

Library  Manager,  DLS-Canberra 

UNIVERSITIES  AND  COLLEGES 

Australian  Defence  Force  Academy 
Library 

Head  of  Aerospace  and  Mechanical  Engineering 
Serials  Section  (M  list),  Deakin  University  Library,  Geelong,  3217 
Hargrave  Library,  Monash  University  (Doc  Data  Sheet  only) 

Librarian,  Flinders  University 


OUTSIDE  AUSTRALIA 

ABSTRACTING  AND  INFORMATION  ORGANISATIONS 

Engineering  Societies  Library,  US 

Documents  Librarian,  The  Center  for  Research  Libraries,  US 

INFORMATION  EXCHANGE  AGREEMENT  PARTNERS 

Acquisitions  Unit,  Science  Reference  and  Information  Service,  UK 

Library  -  Exchange  Desk,  National  Institute  of  Standards  and  Technology,  US 

SPARES  5 


Total  number  of  copies: 


44 


Page  classification:  UNCLASS. 


DEFENCE  SCIENCE  AND  TECHNOLOGY  ORGANISATION 
DOCUMENT  CONTROL  DATA 


FIED 


1.  PRIVACY  MARKING/ CAVEAT  (OF 
DOCUMENT) 


2.  TITLE 

Simulation  of  Passive  Imaging  of  Underwater  Objects 


3.  SECURITY  CLASSIFICATION  (FOR  UNCLASSIFIED  REPORTS 
THAT  ARE  LIMITED  RELEASE  USE  (L)  NEXT  TO  DOCUMENT 
CLASSIFICATION) 


Document 

Title 

Abstract 


(U) 

(U) 

(U) 


4.  AUTHOR(S) 
A.  Kouzoubov 


5.  CORPORATE  AUTHOR 

Aeronautical  and  Maritime  Research  Laboratory 
506  Lorimer  St 

Fishermans  Bend  Victoria  3207  Australia 


6a.  DSTO  NUMBER 
DSTO-TR-1240 


6b.  AR  NUMBER 
AR-01 2-070 


6c.  TYPE  OF  REPORT 
Technical  Report 


7.  DOCUMENT  DATE 
November  2001 


9.  TASK  NUMBER 

10.  TASK  SPONSOR 

11.  NO.  OF  PAGES 

12.  NO.  OF 

NAV99/062 

DNW 

46 

REFERENCES 

35 

8.  FILE  NUMBER 
9505-21-84 


13.  URL  on  the  World  Wide 

http://www.dsto.defence.gov.au/corporate/reports/DSTO-TR-1240.pdf 


14.  RELEASE  AUTHORITY 
Chief,  Maritime  Operations  Division 


15.  SECONDARY  RELEASE  STATEMENT  OF  THIS  DOCUMENT 

Approved  for  public  release 


OVERSEAS  ENQUIRIES  OUTSIDE  STATED  LIMITATIONS  SHOULD  BE  REFERRED  THROUGH  DOCUMENT  EXCHANGE,  PO  BOX  1500,  SALISBURY,  SA  5108 

16.  DELIBERATE  ANNOUNCEMENT 

No  Limitations 

17.  CASUAL  ANNOUNCEMENT _ Yes  . . 

18.  DEFTEST  DESCRIPTORS 

Computerized  simulation.  Imaging  techniques.  Underwater  objects.  Spectral  reflectance.  Air  water 
interactions,  Ray  tracing.  Radiance,  Irradiance,  Mine  hunting 

19.  ABSTRACT 

A  computer  model  simulating  performance  of  passive  optical  sensors  for  imaging  underwater  objects 
from  above-water  paltforms  is  described.  The  model  is  based  on  the  models  of  sun  and  sky  spectral 
radiance,  wavy  sea  surface,  light  propagation  in  seawater  and  is  able  to  render  images  of  underwater 
objects  as  seen  through  the  air-water  interface. 


Page  classification:  UNCLASSIFIED 


